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E(g = CjKE (t ) + C2PE(t n ) for kinetic energy 
and/or potential energy components 
aUa sequence of times {t n } n , where 
C-j and C2 are selected real values, not both 0 


Provide reference values E(t' n ;ref) 
for energy components 
at a sequence of times {t' n } n 


Compute an index comparison 
CHEtg, E(t’ n ;ref) 
for at least one time pair (t n ,tn) 


When the companion index valued lies 
outside a selected range for this index, 
interpret this condition as indicating 
that the estimated kinetic energy 
component and/or potential energy 
component is anomalous 
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Flight Parameter 

Time: t = t n 

FP (t n ;meas) = 

FP (t n ;hist;m)= 

Statistical percentage for FP value = 

Possible source(s) of a typicality = 
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HISTORICAL ANALYSIS OF AIRCRAFT 
FLIGHT PARAMETERS 

CROSS REFERENCE TO RELATED 
APPLICATIONS 

This application is a Continuation In Part of Li.S. Scr. No. 
10/956,523, filed 22 Sep. 2004 now U.S. Pal. No. 7,075,457. 

ORIGIN OF THE INVENTION 

This invention was made, in part, by one or more employ- 
ees of the U.S. government. The U.S. government has the 
right to make, use and/or sell the invention described herein 
without payment of compensation therefor, including but not 
limited to payment of royalties. 

FIELD OF THE INVENTION 

This invention relates to monitoring, analysis and graphic 
illustration of historical energy, location and orientation 
parameters for an aircraft in various phases of a (light. 

BACKGROUND OF THE INVENTION 

An aircraft that is ascending following takeoff or descend- 
ing on approach will have measurable kinetic energy and 
potential energy components, and these components will 
change with time in measurable, if not predictable, manners. 
Desirable energy slates for both takcolT and landing can be 
determined from aircraft manufacturer guidance for these 
phases of flight. For example, where the approach occurs at 
an airport with an operable and reliable instrument landing 
system (ILS), the ILS system may provide data recorded on 
the aircraft to serve as a standard for comparing observed 
kinetic and potential energy components for an aircraft near 
the ground, below 2500 feet altitude and for an assumed 
straight path to a touchdown site. If the airport has no 
operable and reliable ILS, or if the aircraft is not near the 
ground, another mechanism for providing a standard for 
measurements or estimates is needed. On takeoff, where no 
electronic guidance comparable to the glideslope is avail- 
able, the aircraft climb profile can be compared to manu- 
facturer guidance or to observed performance for recorded 
aircraft departures from the particular airport. 

The airline industry has become concerned with the 
problem of unstable aircraft approaches, because approach 
and landing accidents often begin as unstable approaches. 
An “unstable approach” is often defined as an approach 
where below a threshold altitude (1000 feet for IFR and 500 
feel for VFR), the aircraft is not established on a proper glide 
path and with a proper air speed, with a stable descent rate 
and engine power setting, and with a proper landing con- 
figuration (landing gear and flaps extended). Airlines have 
developed approach procedures that call for abandonment of 
an approach that is determined to be unstable. 

Development and testing of methods for delecting atypi- 
cal flights by N.A.S.A. has revealed that high energy during 
an arrival phase (below 10,000 feel but before beginning an 
approach) is the most common reason for a (light to be 
identified as atypical or out of a statistically normal range. 
An atypical high energy arrival phase often corresponds to 
aircraft kinetic energy and/or potential energy that requires 
dissipation of 10-30 percent more energy than is required 
for a reference arrival phase. A reference arrival phase may 
correspond to about a 3 miles per 1000 feet elevation change 
(“3-to- 1 ”) glide path slope and decelerating to an airspeed of 


2 

about 250 knots during descent through 10,000 feet altitude 
to a standard reference speed around 2,500 feet altitude, 
when beginning an approach. 

More than half of the high energy arrivals identified by 
5 atypicality analysis were brought under control within sta- 
bilized approach criteria; some of the remainder of the high 
energy arrivals were abandoned. In contrast, where these 
findings were used to define and search for a high-energy 
arrival exceedance, about three times as many exccdances 
to were detected; and the resulting unstable approaches were 
found to occur more frequently than the recoveries. 

It may be possible to identify, hy historical analysis, a first 
class of high energy arrivals where recovery and subsequent 
stabilization is possible and relatively easy, and a second 
15 class of high energy arrivals in which recovery and subse- 
quent stabilization is likely to be difficult or impossible. 
However, the present procedures for determining presence 
of a reference (acceptable) approach include an electronic 
glide slope that extends linearly from the end of a target 
20 runway to the aircraft, whereas a reference aircraft approach 
path is curved and follows the electronic glide slope only 
from about 1,800 feet above the field to the end of the 
runway. 

A 3-to-l glide path slope, corresponding to decrease of 
25 1 ,000 feet in altitude for every 3 nautical miles horizontal 
travel, is often desirable during an arrival phase. Air speed 
is 250 knots or less by regulation below 10,000 feet, and the 
aircraft decelerates to a lower reference speed before joining 
the approach path. These parameters arc directly available 
30 but are unlikely to prove to be the only relevant parameters 
in determining whether a flight arrival phase is normal or 
other than normal. 

When an energy component value or orientation compo- 
nent value for a completed flight of interest (referred to 
35 herein as a “target flight”) has been measured or observed 
and compared with a corresponding value for a reference 
flight, this information should be displayed for possible 
remedial action on a subsequent flight. A flight operator may 
also benefit from a display of one or more predictions, based 
40 upon the observed or measured target FP values, of the 
behavior of this FP value over a short time interval extend- 
ing into the future. 

What is needed is a system for displaying energy and 
other flight parameters associated with one or more phases 
45 of target flight, which permits historical analysis and visual 
and/or alphanumeric comparison of the target FP values with 
corresponding reference FP values for other flights. Prefer- 
ably, the system should provide corresponding variables for 
a reference flight, for comparison with the target flight, and 
50 should provide a band surrounding of reference FP values 
that indicates values of that FP that are acceptable in 
executing a particular maneuver and ranges of values of that 
FP from which recovery to a reference flight configuration 
is unlikely or substantially impossible. Preferably, a differ- 
55 ence between the target FP value and the reference FP value, 
and one or more time derivatives of this difference should be 
displayed and arc used to predict values of this difference 
over a short time interval in the future. 

60 SUMMARY OF THE INVENTION 

These needs arc met by the invention, which provides a 
method and system for displaying lime variation of one or 
more flight parameter values, including but not limited to 
65 total energy, kinetic energy, potential energy, applied power, 
vertical speed, height above ground, relevant drag indices 
and angle of attack for an aircraft in motion and for variation 
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with time of any of these variables with one or more of 
approximately 20 primary parameters that arise in an energy 
eonfiguration analysis of the aircraft. More particularly, the 
system can compare selected variables for the target flight 
with corresponding variables for a reference (light in a 5 
selected flight phase (c.g., approach to touchdown or take- 
off). Optionally, the system displays target flight parameter 
values and indicates what actions might have been taken 
during the flight to bring the target flight parameter values 
within a percentage band of historical data for the flight to 
parametcr(s). A display of a flight parameter value may be 
graphical, alphanumeric, or a combination of graphical and 
alphanumeric. 

The system displays a percentage band PB including a 
selected percentage value p, in a range such as 15 
70%Sipfi95% of all historic data for a given flight param- 
eter for a similar environment. The system also measures (or 
estimates) and displays a target FP value for a flight of 
interest under similar environmental conditions, for com- 
parison. Optionally, when the target FP value lies outside the 20 
PB the system performs a further analysis to identify what 
anomalies are sources of these conditions. 

In one embodiment, the system measures and analyzes 
relevant parameter values for an ascending or descending 
aircraft to determine if an energy and/or orientation FP value 25 
of the target flight is within, or is outside of, a range for a 
normal flight. This invention can be used in post-flight 
review of flight data and/or as part of a flight operations 
quality assurance program to alert an analyst to presence of 
an anomalous or atypical energy state in historical data. 30 

This measuremcnl/eslimalion/analysis process may 
include the following: 

(i) providing an estimate or measurement of a target flight 
parameter value FP(t„) (referred to as a “measured target FP 
value”) of an aircraft flight parameter during a selected 35 
phase (e.g., takeoff, ascent, descent or approach) of a flight, 

at each of a sequence of measurement (or estimation) times 
t„ (n=l, . . . , N; NS2); 

(ii) providing and displaying a percentage band (“PB”) of 
historical data FP(t„;hist;m) (m=l, . . . , M; Mi£2) for the 40 
flight parameter of interest, drawn from historical FP values 
for M flights under similar conditions; 

(iii) determining if the target FP value for the measure- 
ment time t„ lies within the PB; 

(iv) when the target FP value does not lie within the PB, 45 
visually or aurally indicating this, and optionally recom- 
mending at least one action that may begin to bring subse- 
quently received FP values for the target flight within the PB 
for future measurement times l„; and 

(v) when the target FP value lies within the PB, optionally 50 
displaying this value and the band for the measurement time 

1 ,,- 

Flight parameters that can be monitored, analyzed and/or 
displayed using this approach include: kinetic energy 
KE(t)=m(t)-v(t) 2 /2+coI-co/2; potential energy PE(t)=m(l) 55 
-g--h(t); energy component E(l)=dTKE(t)+d2-PE(t), where 
(dl,d2) are selected non-negative real numbers; energy 
component time derivative (d/dt)E(t), thrust power, vertical 
speed, ground air speed, aircraft mass, height above ground, 
flap position, speed brake position, landing gear position, 60 
other drag indices, roll, pilch and/or yaw angles; and angle 
of attack. 

BRIEF DESCRIPTION OF THE DRAWINGS 

65 

FIGS. 1A and IB illustrate environments in which the 
invention can be practiced. 
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FIGS. 2 and 3 arc flow charts of procedures for analyzing 
aircraft energy components according to the invention. 

FfG. 4 graphically illustrates variation of drag force with 
drag appliance activation. 

FIG. 5 illustrates a display screen that incorporate the 
invention. 

FIG. 6 illustrates variation of a percentage band for 
energy with air miles to touchdown. 

FIGS. 7 and 8 are a How chart and accompanying screen 
for practicing an embodiment of the invention. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIGS. 1A and IB illustrate environments for an ascending 
aircraft (1A) and for a descending aircraft (IB) where the 
invention can be practiced. In FIG. 1A, an aircraft 11A is 
ascending, either after takeoff or in moving from a first flight 
altitude to a second flight altitude. The aircraft has at least 
one of an associated kinetic energy component KE(t„) and/o 
associated potential energy component PE(t„), measured or 
estimated or otherwise provided, at each of a first sequence 
{t„}n of two or more time values, thrust power, vertical 
speed, height above ground, individual or collective drag 
indices, roll angle, pitch angle, yaw angle and angle of 
attack. The aircraft kinetic energy and potential energy 
components for a target flight are, respectively, 

KEU)=m(i)Mifl2+<Mx)-H0Mi)l2, ( 1 ) 

PEU)=m(l)i:hU), ( 2 ) 

where m(t) is the instantaneous aircraft mass (taking account 
of fuel consumption), I(t) is an instantaneous moment of 
inertia tensor for the aircraft, o(l) is an aircraft rotation 
vector, computed with reference to a center of gravity or 
other selected location determined with reference to the 
aircraft (optional), v(t)=dx/dt is the instantaneous aircraft 
velocity and h(l) is the instantaneous height of aircraft eg 
above local reference height, such as local ground height. 
The rotational component of kinetic energy may be negli- 
gible or may be ignored for other rcasonsFor an approach to 
touchdown, the flight parameter of greatest concern is often 
kinetic energy KE(t). 

FIG. 2 is a flow chart of a procedure for practicing an 
embodiment of the invention. In step 21, an aircraft system 
measures or estimates or otherwise provides a value (re- 
ferred to as an “measured value” for convenience herein) 

E(t„)=d) KE(l„}+il2 PE(l„) ( 3 ) 

of an energy component of an aircraft during an ascent phase 
or descent phase of a target flight, at each of a first sequence 
of times (n=l , . . . , N 1 ; N 1 §2), where dl and d2 arc selected 
real values, not both 0. In step 22, the system provides or 
computes a reference value E(l'„;ref) of the energy compo- 
nent at a time, t=t'„, determined with reference to the time t„ 
(n=l, . . . , Nl). The time sequence {t'„} may substantially 
coincide with the sequence {l„}, or each time value t'„ may 
be displaced by a calculable or measurable amount from the 
corresponding time value t„. In step 23, the system computes 
an index of comparison value Cl{E(t„), E(l'„;ref)} of the 
measured and reference energy components for at least one 
time value pair (t„,l'„). When the comparison index value Cl 
lies outside a selected range for this index, the system 
interprets this condition as indicating that the measured 
energy component is anomalous or non-normal or may lead 
to an unstable aircraft maneuver, in step 24. 
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A variety of comparison indices Cl can be used here. 
Some examples arc: (1) a first ratio E(t„)/E(l'„;rcf); 
(2) a second ratio E(i'„;rcf)/E(l„); (3) a difference E(l„)-E 
(t'„;ref)}; (4) an absolute difference IE(l„)-E(t'„;rcf)l; (5) a 
normalized difference {E(l„)— E(l'„;rcr)}/{a-E(t„)+(l-a)-E 
(t'„;rcD}, where a is a selected real value in a range OgaS 1 ; 
(6) a weighted average of the differences KE(t„)-KE(t'„;reO 
and/or PE(t„)-PE(i’„;rcf), such as 


Nl 

IV/t = Yj »'„(! KE(I„) - KE(l'„-, ref)f + I PE(t„) - PE(t’„; re])\’ , \' lp , 


(4) 


where p is a selected positive number (c.g., p=l or 2 or 3.14) 
and {w„}„ is a sequence of weight values (preferably, but not 
necessarily, non-negative); and (7) a monolonic function of 
one or more of the preceding combinations. 

The comparison index Cl may use two or more point 
values, E(l„) and E(t' (1 ;re0, or may use a weighted average 
of these values, such as the average 


Nl 

WM = Y 1 VnCl{E(t«). EiC ref)Y 


(5) 


The analysis may be extended to consider time rates of 
change, (d/dt)KE(l) and/or (d/dl)PE(l), of the kinetic and/or 
potential energy components at a sequence of one or more 
times {t„}„ (n=l, . . . , N2; N211), plus corresponding 
reference lime rates of change, (d/dt)E(t;ref), at a sequence 
of times {t’„} (1 determined with reference to the time 
sequence {l„}„. Another comparison index, C2{ (d/dl)E(t„), 
(d/dt)E(t'„;ref)}, which may be the same or different from 
the comparison index Cl, is computed and compared with a 
second selected range to determine if the aircraft flight is 
anomalous or non-normal or is within a normal range. 
Again, the comparison index C2 may use point values or a 
weighted average of the values (d/dt)E(l) and/or (d/dt)E(t; 
ref). 

FIG. 3 is a (low chart of another procedure for practicing 
an embodiment using time derivatives of the energy com- 
ponent E(l„). In step 31, an aircraft system measures or 
estimates or otherwise provides an “estimated value”, (d/dt) 
E(t„)=d3 (d/dt)KE(t„)+d4.(d/dt)PE(t„) of an energy compo- 
nent of an aircraft during an ascent phase or descent phase 
of a target flight, at each of a first sequence of times 
(n=l, . . . , Nl; Nl §2), where d3 and d4 are selected real 
values, not both 0. In step 32, the system provides or 
computes a reference value (d/dl)E(l'„;ref) of the energy 
component at a time, t=t'„, determined with reference to the 
time t„ (n=l, . . . , Nl). The time sequence {t"„} may 
substantially coincide with the sequence {t,,} or {t'„}„, or 
each time value l"„ may be displaced by a calculable or 
measurable amount from the corresponding time value t„. In 
step 33, the system computes an index of comparison value 
C2{(d/dt)E(t„), (d/dl)E(l'„;rcf)} of the estimated and refer- 
ence energy component time derivatives for at least one time 
value pair (t„,l’’„). When the comparison index value C2 lies 
outside a selected range for this index, the system interprets 
this condition as indicating that the estimated energy com- 
ponent time derivative is anomalous or non-normal or would 
lead to an unstable aircraft maneuver, in step 34. 

The analysis may be further extended to consider a third 
comparison index, C3{E(t H ), E(l'„;rcf), (d/dt)E(l H ), (d/dt)E 


(l"„;rcf)}, that depends upon some or all of the estimated 
values and time rates of change of the estimated values of 
the energy components. Again, the comparison index C3 
may use point values or a weighted average of the values 
5 E(t) and/or E(t;ref) and/or (d/dt)E(t)/di and/or (d/di)E(l;rcO. 

A formulation of, and use of, the equations of motion of 
a target aircraft flight, including the effects of gravity, 
variable wind speeds, drag and lift forces on various control 
surfaces, variation of aircraft mass due to fuel consumption, 
to and variable thrust, is set forth in an Appendix. A thrust 
vector is determined, as a function of the location coordi- 
nates, that will move the aircraft from an initial velocity 
vector v H (x 0 ,y 0 ,z 0 ) to a desired final velocity vector v/x^y^ 
z r ) as part of a takeoff phase or as part of an approach phase 
15 for a flight. The aircraft kinetic energy is a sum as in Eq. (1). 

Each aircraft has an associated group of drag indices, one 
for each activalablc drag appliance (landing gear, wing llap, 
spoilcr/specd brake, etc.). Each drag index has a maximum 
value where the drag appliance is fully activated and has a 
20 spectrum of drag values extending from zero activation 
through partial activation to full activation of the appliance, 
as illustrated schematically in FIG. 4. With the drag appli- 
ance completely inactivated, the corresponding drag index is 
normally 0. The drag force associated with one drag appli- 
25 ancc is assumed to be independent of the drag force asso- 
ciated with another drag appliance, in a first approximation. 
In an approach to landing, for example, where a relatively 
small amount of additional drag force is required for fine 
adjustment, one or more drag appliances can be partly or 
30 fully activated to provide this additional drag force, relying 
on information illustrated in FIG. 4 for each drag appliance. 
If the amount of additional drag force needed for the 
adjustment is greater than the maximum drag force associ- 
ated with all the activated drag appliances, the aircraft would 
35 need to use additional procedures to provide the additional 
drag force, or the target flight configuration should be (or 
should have been) terminated and reconfigured. In practice, 
some drag appliances, such as landing gear, are normally 
inactivated or fully activated, while other drag appliances, 
40 such as a speed brake, have a near-continuous range of 
settings. The sum of the drag indices for all (activated) drag 
appliances is determined and provided s a supplement to the 
drag force(s) provided by the other aircraft components. 

Monitoring of thrust power developed by the engine(s) of 
45 the aircraft is straightforward and is an important control 
variable in change of the energy component E(l) defined in 
Eq. (3). Thrust developed can be estimated using measured 
fuel flow rate, temperature within the engine(s) and other 
relevant variables. 

50 Aircraft angle of attack can be measured, made available 
and recorded on the aircraft. 

FIG. 5 is a characterization of a sequence of measured 
values for a flight parameter FP(l) of interest over a selected 
time interval, for the target flight, and for historical data 
55 from a collection of preceding flights under similar condi- 
tions, denoted FP(t;hist;m) for M preceding flights 
(m=l, . . . , M; MS2). The time l is referenced to a sequence 
of measurement limes t„, corresponding to a reference time, 
such as the time before lift-off, the time since lift-off, the 
60 time since passing a specified waypoint, or the time preced- 
ing touchdown of the aircraft. For any measurement time t„, 
the historical data FP(t„;hist;m) provide a spectrum of FP 
values, as illustrated in FIG. 5, which can be arranged from 
lowest to highest as shown. A percentage band PB of 
65 consecutive FP values (from lowest to highest, or from 
highest to lowest) for a corresponding measurement time t„ 
is specified, corresponding to a selected percentage in a 
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range, such as 70%SpS95%, and the collection of FP 
values in this band is used as a standard against which the 
target FP value FPd) is compared. For example, the per- 
centage band PB may be all values FP(ln;hisl;m) in (1) the 
lowest p percent, (2) the highest p percent; (3) a symmetric 
band, centered at the median value; or (4) a band having the 
lowest max-min difference, FP(max)-FP(min), for all val- 
ues in the band. 

The flight parameter value FP(t„) for the target flight is 
then compared with values FP(t„;hisl;m) in the band PB, as 
illustrated in FIG. 5, to determine if the target FP configu- 
ration is within a normal or typical range of the historical FP 
data for the time l=t„. Where the target FP value lies outside 
the percentage band PB, this target FP value is optionally 
interpreted as atypical or non-normal, and a recommenda- 
tion for specified action is optionally provided to bring this 
target FP value within a typical range. A target FP value may 
lie in a normal range for a first sequence of measurement 
times l„ and may lie in an atypical range in a second 
sequence of measurement times t„. The PB may be charac- 
terized by upper and lower traces representing maximum 
and minimum FP values within the band, as shown in FIG, 
5 , or another FP display scheme may be used. Whatever 
characterization of the PB is used, the target FP values FP(t„) 
for each relevant measurement time t„ are explicitly dis- 
played on the same graph for monitoring by the aircraft 
operator. 

Where the target FP value FP(t„) lies substantially outside 
the percentage band PB for that measurement lime, the 
system optionally performs a further analysis to (i) indicate 
presence of an atypical or anomalous FP value; (ii) estimate 
a percentage band (e.g., highest or lowest 1 .5 percent in the 
statistical polulalion of values for that FP) in which the FP 
value falls; and/or (iii) identify one or more sources of the 
anomalous value. 

For example, where the kinetic energy of the target 
aircraft on approach was too large relative to the FP values 
within the PB, indicating that the approach velocity was too 
large, the system may identify a kinetic energy value for a 
preceding waypoint or for a preceding altitude during 
descent that was much higher than an acceptable value. 
These recommendations arc made for an already-completed 
flight but may be useful in comparing similarly atypical 
target flights that are completed. 

FIG. 6 illustrates variation with air miles to touchdown of 
a percentage band and median value for an energy index for 
a normal or reference flight. Note that the PB width 
decreases steadily as touchdown is approached. 

FIG. 7 is a flow chart of a procedure for practicing the 
invention. In step 71, the system receives or otherwise 
provides a target FP value, measured or otherwise provided, 
for a measurement time l„. In step 72, the system provides 
and optionally displays historical data FP(t„;hist;m) 
(m=l, . . . , M) for M preceding flights under similar 
conditions for the FP of interest, displays the target FP value 
and the percentage band PB for the historical FP data 
corresponding to the time t„. Optionally, the system displays 
one or more target FP values received (at measurement times 
t„, <t„) before the target FP value was received and the 
corresponding PB for these previously received FP values. 
In step 73, the system determines if the target FP value is 
within the PB for the corresponding time t„. If the answer to 
the query in step 73 is “yes,” the system takes no further 
action and returns to step 71. If the answer to the query in 
step 73 is “no,” the system indicates, in step 74, that the 
target FP value is atypical, displays FP(t„) and FP(l„;hist;m), 
returns to step 71, and optionally recommends at least one 


8 

corrective action, if any, that could have been taken to bring 
the subsequently received target FP values within the PB for 
at least one future measurement time. 

FIG. 8 illustrates a screen that may be used in connection 
5 with step 74 of FIG. 7, displaying one or more of: the flight 
parameter, the measurement time t„, the presently measured 
FP value FP(t„;mcas), a corresponding historical (ideal or 
reference) FP value FP(t„;hist;m) for a preceding (light 
number m (=1, . . . , M) under similar conditions, an 
to estimated statistical percentage associated with the atypical 
FP value, a possible sourcc(s) of the atypicality, and a 
recommended corrective action, if any (optional). 

What is claimed is: 

1. A method of displaying flight parameters of an aircraft 
15 in flight, the method comprising: 

providing an estimate or measurement of a value (referred 
to as an “estimated value”) of a present flight parameter 
value, FP(t„) for an aircraft during a selected phase of 
a flight, at each of a first sequence of times, t=t„ 
20 (n=l, . . . , N; Ng2), 

providing and displaying a percentage band of historical 
data FP(t„;hist;m) (m=l, . . . , M; Mg2) for the flight 
parameter of interest, drawn from FP values for M 
flights under similar conditions and including no more 
25 than a selected percentage p of the M flights, with 
p<100 percent; 

determining if the present FP value for the measurement 
lime t„ lies within the percentage band; and 

when the present FP value does not lie within the per- 
30 cenlage band, displaying the present FP value and the 
percentage band for the measurement lime t„ in at least 
one of a visually perceptible and an audibly perceptible 
manner. 

2. The method of claim 1, further comprising: 

35 when said present FP value docs not lie within said 
percentage band, taking at least one of the following 
actions: (i) acknowledging that said present FP value is 
atypical or anomalous; (ii) providing an estimate of a 
statistical percentage, among a population of values, 
40 corresponding to said present FP value; and (iii) iden- 
tifying at least one source of said atypical present FP 
value. 

3. The method of claim 1, further comprising: 

when said present FP value lies within said percentage 
45 band displaying said present value and said percentage 
band for said measurement time t„. 

4. The method of claim 1, further comprising choosing 
said percentage band to include a consecutively, monoioni- 
cally increasing group of values FP(t„;hist;m) for said lime 

50 t„. 

5 . The method of claim 4 , further comprising choosing 
said selected percentage p in a range 70 percent<p<95 
percent. 

6. The method of claim 1, further comprising selecting 
55 said percentage p from the following group: (1) the lowest 

p percent, (2) the highest p percent; (3) a symmetric band, 
centered at the median value; or (4) a band having the lowest 
max-min difference, FP(max)-FP(min), for all values in the 
band. 

60 7. The method of claim 1, further comprising choosing 

said flight parameter value from the group of flight param- 
eter values consisting of: an energy component, 

E(t, 1 )=dl KE(t„)+d2-PE(t„), of a combination of a kinetic 
energy component KE(t„) and a potential energy component 
65 PE(t (1 ) of said aircraft, where dl and d2 arc selected real 
values, not both 0; an estimated value (d/dl)E(t„)=d3-(d/dt) 
KE(t, I )+d4-(d/dt)PE(t, 1 ) of a time rate of change of said 
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estimated energy component, where d3 and d4 arc selected 
real values, not both 0 ; thrust power applied to said aircraft; 
weight of said aircraft; flap position for said aircraft; speed 
brake position for said aircraft; a selected drag index for at 
least one drag appliance for said aircraft; air speed ol said 
aircraft; vertical speed of said aircraft; height of said aircralt 
above a local reference height; roll angle of said aircraft; 
pitch angle of said aircraft; yaw angle of sad aircraft; and 
angle of attack of said aircraft. 

8 . A system of displaying flight parameters of an aircraft 
in flight, the system comprising a computer that is pro- 
grammed: 

to provide an estimate or measurement of a value (re- 
ferred to as an “estimated value”) of a present flight 
parameter value, FP(l„) for an aircraft during a selected 
phase of a flight, at each of a first sequence of times, 
t=t B (n=l, . . . , N; N£2), 

to provide and display a percentage band of historical data 
FP(t ( 1 ;hist;m) (m=l, . . . , M; Mi?2) for the flight 
parameter of interest, drawn from FP values for M 
flights under similar conditions and including no more 
than a selected percentage p of the M flights, with 
pel 00 percent; 

to determine if the present FP value for the measurement 
time t„ lies within the percentage band; and 

when the present FP value does not lie within the per- 
centage band, to display the present FP value and the 
percentage band for the measurement lime t„ in at least 
one of a visually perceptible and an audibly perceptible 
manner. 

9. The system of claim 8 , wherein said computer is further 
programmed so that: 

when said present FP value does not lie within said 
percentage band, performinmg at least one of the 
following actions: (i) acknowledging that said present 
FP value is atypical or anomalous; (ii) providing an 
estimate of a statistical percentage, among a population 
of values, corresponding to said present FP value; and 
(iii) identifying at least one source of said atypical 
present FP value. 
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10. The system of claim 8 , wherein said computer is 
further programmed so that: 

when said present FP value lies within said percentage 
band to display said present value and said percentage 

5 band for said measurement time t„. 

11. The system of claim 8 , wherein said computer is 
further programmed to choose said percentage band to 
include a selected percentage of said historical data, where 
the selected percentage includes a consecutively, monotoni- 

10 cally increasing group of values FP(t„;hist;m) for said time 

U 

12. The system of claim 8 , wherein said computer is 
further programmed to choose said selected p in a range 70 
pcrccnl<p<95 percent. 

13. The system of claim 8 , wherein said computer is 
further programmed to select said percentage p from the 
following group: ( 1 ) the lowest p percent, ( 2 ) the highest p 
percent; (3) a symmetric band, centered at the median value; 

20 or ( 4 ) a band having the lowest max-min difference, 
FP(max)-FP(min), for all values in the band. 

14. The system of claim 8 , wherein said computer is 
further programmed to choose said flight parameter value 
from the group of flight parameter values consisting of: an 

25 energy component, E(t, 1 )=dl-KE(l„)+d2-PE(t„), of a combi- 
nation of a kinetic energy component KE(t„) and a potential 
energy component PE(t„) of said aircraft, where dl and d2 
are selected real values, not both 0 ; an estimated value 
(d/dt)E(t„)=d3-(d/dt)KE(t, 1 )+d4-(d/dt)PE(t, I ) of a lime rate of 

30 change of said estimated energy component, where d3 and 
d 4 arc selected real values, not both 0 ; thrust power applied 
to said aircraft; weight of said aircraft; flap position for said 
aircraft; speed brake position for said aircraft; a selected 
drag index for at least one drag appliance for said aircraft; 

35 air speed of said aircraft; vertical speed of said aircraft; 
height of said aircraft above a local reference height; roll 
angle of said aircraft; pitch angle of said aircraft; yaw angle 
of sad aircraft; and angle of attack of said aircraft. 

* * * * * 



